The effect of dynamical Bloch oscillations on optical-field-induced current in a wide 

gap dielectric 
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We consider the motion of charge carriers in a bulk wide gap dielectric interacting with a few-cycle 
laser pulse. A semiclassical model based on Bloch equations is applied to describe the emerging time 
dependent macroscopic currents for laser intensities approaching the damage threshold. At such laser 
intensities, electrons can reach edges of the first Brillouin zone even for electron-phonon scattering 
rates as high as those known for Si02. We find that, whenever this happens. Brag-like reflections of 
electron waves, also known as Bloch oscillations, affect the dependence of the charge displaced by 
the laser pulse on its carrier-envelope phase. 



I. INTRODUCTION 



The motion of conduction electrons in a crystalline 
solid is usually considered to be similar to that in free 
space, apart from scattering processes. The situation is 
radically different in the case vifhere an external electric 
field is so strong that, in spite of scattering, an electron 
can acquire such a high quasimomentum that it reaches 
an edge of the first Brillouin zone. While the kinetic en- 
ergy of a free electron exposed to a constant external field 
would indefinitely increase, an electron in a crystal first 
slows down until it reaches the top of the energy band, 
and then it moves in the opposite direction towards the 
bottom of the band. In the semiclassical picture neglect- 
ing scattering, the electron would move periodically back 
and forth between its initial and final positions. This 
phenomenon is known as Bloch oscillations [1], and it 
leads to Wannier-Stark localization [^] . In the reciprocal 
space, an electron reaching an edge of the first Brillouin 
zone continues its motion from the opposite side of the 
zone. In the real space, this corresponds to a Bragg-like 
reflection of an electron wave jSj] • While Bloch oscillations 
and Wannier-Stark localization are usually considered in 
the case of a constant external field, essentially the same 
physical phenomena take place if the external field is time 
dependent. In this paper, we use the term 'dynamical 
Bloch oscillations' (DBOs) to describe phenomena that 
occur whenever an electron wave is refiected at an edge 
of the Brillouin zone. 

Until recently, Bloch oscillations were thought to be 
impossible to observe in bulk solids because of scatter- 
ing. The period of Bloch oscillations in a constant field 
F is given by Tb = h{eFa)~^, where h is the Planck con- 
stant and a is a lattice period. To observe Bloch oscilla- 
tions, Tb must be smaller than characteristic scattering 
and dephasing times, which are usually on the order of 



Ts ~ 10 ^■^ s. This implies that the external field must 
be stronger that F > h{eaTs)~^ ~ 10® Vm^^. Such 
a strong constant field would destroy even wide gap di- 
electrics. Therefore, over the last few decades, Bloch os- 
cillations were pred ominantly studied in semiconductor 
super lattices |4|-[l0| and optical lattices 11 1. 
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The situation has recently changed as intense few-cycle 
pulses were generated in the mid-infrared (MIR) spectral 
region [l3 , and the duration of the shortest near-infrared 
(NIR) pulses approached one optical cycle [T§|. Ghimire 
et al. [3, fl5| observed that anharmonicity in the motion 
of charge carriers created and driven by an intense MIR 
pulse in ZnO resulted in the generation of high-order har- 
monics [l3 and a red-shift of absorption edge. In their 
parameter regime, the field was intense enough to drive 
conduction electrons beyond the first Brillouin zone, so 
that Bloch oscillations were suggested to be responsible 
for the observed effects. Very recently, Schiffrin et al. 
[13] have found that a 4 fs NIR pulse with a peak elec- 
tric field of 20 GVm~^ can induce measurable currents 
in a Si02 sample. Furthermore, it has been found that 
these currents can be steered by controlling the carrier- 
envelope phase (CEP) [l^ of laser pulses. These findings 
were interpreted in terms of Wannier-Stark states [l^ . 

While it remains debatable whether Bloch oscillations 
played a major role in these particular measurements, 
there is no doubt that intense few-cycle pulses enable 
experiments in the parameter regime where electrons in 
a bulk solid are accelerated beyond the first Brillouin 
zone. The purpose of this paper is to study some basic 
effects related to this parameter regime. In particular, we 
consider the role of electron scattering and dephasing. 

Electron-phonon scattering rates are known to be par- 
ticularly high for Si02 due to a strong coupling be- 
tween conduction electrons and longitudinal optical (LO) 
phonons [13, [HI • At the same time, scattering rates for 
moderately hot (i^kin ~ 1 eV) electrons, which we con- 
sider in this paper, have never been measured directly, 
and there are still open questions related to the scatter- 
ing of very hot conduction electrons [23 ] . 
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II. SYSTEM AND MODEL 

We consider the following model (see Fig. [Ij: the xy 
{z = 0) plane is the surface of a dielectric. Short pulses 
with a stabilized CEP propagate along the z-axis and 
impinge on this surface. In our calculations, we as- 
sume that the external electric field within the sample 
is linearly polarized and given by £x{t) = cos(ajot + 
</5cEp) exp(-^), Ey^ Ez = 0. 




Figure 1: A schematic view of the creation and drivin g of 
macroscopic currents with a laser pulse. Measurements [131 
yield the charge transferred by the pulse. 



We adopt the two-band approximation and consider 
the electron motion in two spatial dimensions. At each 
moment t, we describe electronic excitations in the sam- 
ple with the aid of quantum-mechanical density matrices 



p(k,t) 



nc(k,t) P(k,t) 
P*(k,t) nv(k,t) 



(1) 



where ric, correspond to the conduction and valence 
band populations, the off-diagonal element P(k, £) repre- 
sents the interband coherence, and the quasimomentum 
k has two components: k = (kx,ky). The time depen- 
dence of the density matrix is formally described by 



9 , f dp 



dp 

dt 



dp 

field \ / scatt 



(2) 



describe photoexcitation, the acceleration of charge car- 
riers by the external field and electron scattering. 

We chose to describe photoexcitation dipole transitions 
using the following form of k resolved optical Bloch equa- 
tions Hill: 



ap(k) 

dt 

ducCk) 

dt 
dn^ik) 

dt 



-[Sc(k)-K(k)-iH^(k) 
n 

iK(k) -nv(k)]4v(k)£:(t), 

23[dev(k)£(i)P*(k)], 

dn,{k) 
dt ' 



(3) 



where £'v(k) and £'c(k) are the energies of the valence 
and conduction bands, respectively. 

The phenomenological rate n describes the decay of 
the interband coherences. For simplicity, we estimate the 
dependence of the dipole matrix elements on k as [2^ 



dcv(k) =dcv(0) 



The actual value of c?cv(0) has little qualitative effect on 
our results, as long as saturation-related phenomena are 
negligible, i.e. the excited population is well below unity. 
In the following, we use dcv(O) = 0.1 atomic units. Note 
that Eqs. ([3]) make no use of the rotating wave approxi- 
mation. This allows us to investigate dynamics that un- 
fold within a single optical oscillation of the laser pulse. 
(See e.g. [25| for a detailed discussion of related phenom- 



ena.) 

The external electric field not only causes transitions 
between the bands, but it also accelerates and deceler- 
ates charge carriers. These field-driven dynamics are ac- 
counted for by the second term in Eq. We neglect 

off-diagonal terms in ( ^ ) and evaluate the diagonal 
ones as ^ ^^'^''^ 



/ dn 



= (k) 



dt 



field 



= — -f (t)Vknv,c(k). 

h 



(4) 



The three terms on the right-hand side of this equation 



Note that £{t) includes the screening field due to the 
collective electron response [l^ (i.e. it denotes the field 
that actually interacts with individual electrons). 

The third term in Eq. ^ accounts for the loss of in- 
traband coherence due to scattering, where longitudinal 
optical (LO) phonons are considered to play the major 
role. The electron-phonon interaction is described by the 
Frohlich Hamiltonian, and standard methods [l^l lead to 
the following dynamical equations: 
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'9nc(k) 
dt 



7o 



- 70 ^ 6{Eik + q)-^c(k)-fiwLo) 4 \ -^q"c(k) [1 - ne(k + q)] + (7Vq+l)ne(k + q) [1 - ne(k)] j> + 

catt q ^ 

^5(i?(k - q) - i?e(k) + fi^Lo)^|-(A^q 

n ^ L 



+ l)n,(k) [1 - 7ic(k - q)] + iVq7ic(k - q) [1 - nc(k)] \, (5) 



where ftwLO denotes the energy of an LO phonon, which 
is assumed to be independent of the reciprocal-space vec- 
tor q. The related phonon density is denoted by N^, and 
7o is an electron-phonon coupling constant. For the sake 
of simplicity, we use a tight-binding type dispersion rela- 
tion, i.e. i?c(k) is proportional to 2 — cos(fca;a) — cos(fcj,a), 
where a is a lattice constant. At room temperature, the 
average thermal phonon number and, consequently, A'q 
are practically zero, so phonon emission processes (ac- 
companied by electron scattering events with a loss of 
electron energy) dominate the scattering dynamics. 

The model above allows us to calculate the rate of scat- 
tering on LO phonons 7(k), which is proportional to 70 
appearing in Eq. ([5]). This scattering rate determines 
how fast the quasimomentum of an electron wave packet 
decreases as a consequence of scattering events. Our nu- 
merical calculations show - in accordance with the ana- 
lytical results presented e.g. in ^ - that the rate 7(k) 
has a pronounced minimum at |fc| = 0, being nearly con- 
stant (7(k) « 7) for kinetic energies in the range between 
0.5 eV and 2 eV. In the following, we use 7 as a label 
to quantify the strength of electron-phonon interaction 
in different simulations, although we used k dependent 
scattering rates in our calculations. 

For our simulations, we used material parameters that 
correspond Si02 : a band gap of 9 eV, lattice period a = 
0.5 nm, two LO phonon modes with energies ftui^o = 
0.153 and 0.063 eV, and a combined scattering rate equal 
to 7 = 0.3 fs^^. The laser pulse parameters correspond 
to the experiment described in [l7|: ujq — 2.51 fs~^ and 
r = 2.3 fs. 



III. RESULTS 

A. Electron acceleration in the laser field 

Before we present simulations where a laser pulse cre- 
ates and drives charge carriers, let us consider the laser- 
driven motion of an initially free electron neglecting in- 
terband transitions. Similar simulations can be found 
in As an example, we take the initial distribution 

of conduction electrons as a Gaussian wave packet cen- 
tred at k = 0, neglect the terms (§f ) and (§f ) 

V / cxc V / scatt 

in Eq. ([2]) and model the time evolution of the electron 
wave packet by solving Eq. ^ . We plot the distribution 
of conduction electrons in false-colour diagrams, where 
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Figure 2: Time evolution of a wave packet initially centred 
at Ic = in a simulation that neglects photoexcitation and 
electron-phonon scattering, (a) The distribution of conduc- 
tion electrons nc(k,t) in the ky = plane. The dashed white 
line is the 'semiclassical trajectory' evaluated with the aid 
of the acceleration theorem ([6]). (b) The electric field £x{t) 
that drives the wave packet, the current density jx{t) and the 
time dependent transferred charge Q{t). The parameters are 
£0 = 4 GVm-\ <^CEP = 0, T = 2.3 fs. 



colours vary from black through red to yellow as the 
electron density nc(k, i) grows from zero to its maximal 
value. The top panel of Fig. [5] presents such a diagram 
in the plane ky = 0. In this example, the distribution 
remains localized in the reciprocal space, and it is dy- 
namically shifted by the field of the laser pulse. The 
reciprocal-space motion of the wave packet is appropri- 
ately described by the 'acceleration theorem' [2^ : 



'dt 



-Sit). 



(6) 



A solution of this equation with the initial condition 
fc(imin) = is shown by the dashed white line in Fig. [21 
This solution can be regarded as a trajectory of a 'classi- 
cal particle' in the reciprocal space. Let us note that as 
long as neither scattering nor excitation is taken into ac- 
count, there is a simple scaling in the model: increasing 
both the carrier frequency ijJq and the amplitude of the 
laser pulse £q by a certain factor is equivalent to choosing 
a new unit of time in Eq. and it does not change the 
maximal quasimomentum that a wave packet can reach 
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Figure 3: Time evolution of a wave packet initially centred 
at k = in a simulation that neglects photoexcitation but 
takes electron-phonon scattering into account. 7 = 0.3 fs~^, 
all other parameters are the same as in Fig. (2] 



in the reciprocal space. 

Having evaluated the time evolution of the density ma- 
trix, we are to investigate measurable physical quantities. 
The time dependent current density j{t) is a sum of con- 
tributions from conduction-band electrons 

jcW = / evc(k)ne(k,0d'k, (7) 
Jbz 

valence-band holes 

jv(t)= / evv(k)nv(k,t)d2k, (8) 
Jbz 

and interband coherences 

jcv(i) f P(k,t)4v(k) [S„(k) - E,{k)] d^k, (9) 

n Jbz 

where the integrals are taken over the first Brillouin zone, 
and the velocity distributions are given by Vv,c(k) — 
|-Vki?v,c(k). Eqs. ([T])-® result from evaluating the ex- 
pectation value of the current operator averaged over 
a unit cell for a state described by the density matrix 
p^. Integrating the current density with respect to time, 
we obtain the charge that flows per a unit-area surface 
perpendicular to the x-axis: Q{t) — jx{t') dt' = 
Qcit) + Qv{t) + Qcv{t) , where jj; denotes the x component 
of the total current density j = jc+jv+jcv It is the charge 
displaced by the laser pulse Q = Q{oo) that can be mea- 
sured in experiments [17|. In the following, we assume 
that jc gives the dominant contribution to Q. Indeed, jv 
is negligible in comparison with jc due to the low mobility 
of holes in the valence band, while jcv mainly describes 
the polarization response of valence-band electrons; ac- 
cording to our calculations, Qcv{t) is proportional to the 
applied field within a relative error of no more than 5% 



up to a field amplitude of £0 = 25 GVm~^, Qcv(oo) be- 
ing negligibly small in comparison to Qc{oo). With this 
in mind, we restrict our analysis to jc and Qc, referring to 
them as 'current density' and 'transferred charge', respec- 
tively. These quantities are shown in Fig. [2jb), together 
with the electric field of the laser pulse, which is depicted 
by the dashed black line. 

In Fig. 131 we show the effects of electron-phonon scat- 
tering on the dynamics of a conduction-band wave packet 
initially centred at k = 0. A comparison with Fig. [5J 
where scattering was neglected, reveals that the electron- 
phonon interaction slightly reduces jx{t)- This is ex- 
plained by the fact that electrons get scattered out of 
the plane ky = 0. 

(a) I , , , , , , , 1 
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Figure 4: Average conduction-band population (ric) = 
(/d^k) J nc(k)d^k for different interband coherence re- 
laxation rates k. For k = 5.0 fs~^ (dashed blue curve), 
we plot 0.05(nc). The parameters are £0 ~ 5 GVm~^, 
fCEP = 0,r = 2.3 fs, 7 = 0.3 fs'^ 



B. Nonresonant interband excitations 

In order to expose a medium to a very intense laser 
field without destroying it, the medium must be possi- 
bly transparent, so that little energy will remain in the 
medium after the interaction with the pulse. This implies 
that the laser frequency wq should be much smaller than 
the bandgap. This is why we are interested in studying 
nonresonant excitations. With our laser parameters, it 
takes more than five laser photons to bridge the bandgap 
of SiOa. 

The dynamics of interband excitations predicted by 
our model are shown in Fig. 2] For k — 0, Eq. ([3]) de- 
scribes a completely coherent excitation process, with the 
average population in the conduction band being roughly 
proportional to the laser intensity. These periodic exci- 
tations and deexcitations are sometimes referred to as 
'virtual excitations', as they largely represent distortions 
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of initial electronic states. Both virtual and real excita- 
tions contribute to photocurrents [l^. In the opposite 
extreme, if we assume an unrealistically fast decoherence 
K = 5 fs~^, the conduction band population becomes a 
monotonically increasing step-like function of time. Lit- 
tle is known about ultrafast dephasing in Si02, so we 
use K = 0.1 fs"^ for our further simulations as a value 
that corresponds to an intermediate regime of interband 
excitations. 



C. Laser-driven motion of photoexcited charge 
carriers and the effect of the carrier— envelope phase 

Here, we investigate outcomes of our model with all 
three terms of Eq. ^ being taken into account. The 
most interesting outcome of combining excitation and 
laser-driven motion is that this may result in a nonzero 
charge Q transferred by a laser pulse — a transport effect 
that does not occur in simulations neglecting interband 
transitions, as Figs. [2] and |3] illustrate. In this section, 
we show that this is a combination of interband transi- 
tions and a dispersion law that determines the transferred 
charge. 




-4 -2 2 4 6 

Time {fs) 

Figure 5: (a) Time evolution of the conduction-band elec- 
tron population that emerges due to multiphoton excitations. 
The cross section in the kx direction (parallel to the ex- 
ternal field polarization) is shown for the following param- 
eters: f = 3 GVm~^,(pcEP = 0,r = 2.3 fs, k — 0.1 fs^^, 
7 = 0.3 fs~^. The dashed white line represent a semiclassical 
trajectory released at i = with k = 0. (b) The electric field 
of the laser pulse £x{t), the induced current density jx{t) and 
the transferred charge Q{t) as functions of time. 

In the upper panels of Figs. [5] and [51 we show nc{kx,t) 
obtained in simulations that account for all the relevant 
processes: multiphoton excitations, light-driven motion 
of charge carriers, dephasing and electron-phonon scat- 
tering. The simulation were performed for a 'cosine' 
(VCEP = 0) and a 'sine' (vscep = '"'/S) laser pulses, 
respectively. The lower panels of the figures show the 



electric field of the laser pulse £x {t) , the induced current 
density jx{t) and the transferred charge Q{t) as functions 
of time. Because of scattering, the current density at the 
end of each simulation quickly approaches zero, but the 
final value of the transferred charge is in general nonzero, 
and it has a significantly higher value for the cosine pulse 
(Fig.©. 

The dashed lines in Figs.[Sl[a) andlHIa) represent 'semi- 
classical electron trajectories' released at a peak of the 
electric field, which are solutions of Eq. ^ with the ini- 
tial condition k(io) — 0, the initial time being to = for 
the cosine pulse and to — ~tt{2ujo)^^ for the sine pulse. 
Even without taking scattering into account, this semi- 
classical analysis can be used to explain many features 
observed in our calculations. In this picture, the con- 
tribution from a particular electron to the transferred 
charge at a final time ^max is determined by the semi- 
classical electron displacement: 

s{to)= f '^^\,{k{t)) dt, (10) 

where to is a time when the electron appeared in the con- 
duction band, k(t) satisfies Eq. Vc(k) = ■iVk-Ec(k) 
and an implicit assumption was made that the initial 
velocity of the electron is zero. In the case of a short 
cosine pulse, the central half-cycle of the electric field 
has a significantly higher amplitude than any other half- 
cycle. Consequently, there is one dominant semiclassi- 
cal trajectory that starts at the peak of the main half- 
cycle. A simple calculation shows that, in our example, 
the semiclassical displacement associated with this tra- 
jectory is negative, which agrees with the positive final 
transferred charge in Fig. [SJ For a short sine pulse, there 
are two dominant trajectories which start at the peaks of 
the two most pronounced half-cycles of the electric field. 
The electron displacements associated with these trajec- 
tories have opposite signs and 'add destructively' in the 
case shown in Fig. [B) This explains why the magnitude 
of the net transferred charge is considerably smaller for 
<y5cEP = 7r/2 than for ipcEP — 0. 

We further elaborate on the role of the CEP in Fig. [71 
were we plot Q{(pcEp) for different values of k and 7. 
From this figure, one can see that (5(v?cep + tt) — 
— (3((/3CEp), which is a direct consequence of symmetry: 
adding tt to the CEP is equivalent to substituting £{t) 
with —£{t), which is equivalent to replacing x with —x in 
our symmetric arrangement. One also infers from Fig. [7] 
that the relaxation rate of interband coherences infiu- 
ences the positions of minima and maxima of Qi^pcEp), 
which is not the case for the phonon scattering rate 7. 

Increasing the amplitude of the input laser fields, we 
reach the regime where even electrons initially excited in 
the middle of the Brillouin zone experience refiections at 
its edges during the laser pulse — dynamical Bloch oscil- 
lations (DBOs) take place. This is illustrated in Fig. |8l 
where we plot the outcomes of a simulation with a laser 
pulse that has a peak amplitude of the electric field of 
10 GVm~^, the other parameters being the same as in 
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Figure 6: The same as Fig. [5l but for a sine pulse (</5cep ~ 
it/2). The semiclassical trajectory begins at the first main 
maximum of the laser field intensity, i.e. at t = — 7r(2a;o)^^ ~ 
0.6 fs, and it trespasses k = at the peak of the next half- 
cycle, which is the starting point of the second dominant elec- 
tron trajectory. 




Figure 7: The dependence of the total charge transferred by a 
laser pulse on its carrier-envelope phase (CEP). The param- 
eters are £o = 2 GVm~^, r — 2.3 fs, 7 = 0.3 fs~^. The curve 
corresponding to k = 0.1, 7 = 0.05 fs~^ has been normalized; 
the other ones have been scaled by the same number, i.e. the 
amplitudes of the curves relative to each other are kept fixed. 
Note that we use k = 0.1, 7 = 0.3 fs~^ in the following. 



the previous simulations. DBOs occur in spite of the 
fact that we use a fairly large value for the the electron- 
phonon scattering rate (7 — 0.3 fs~^). In contrast to the 
case of a low field (Fig. [SJ , the final transferred charge is 
now negative, which can be interpreted in terms of the 
semiclassical electron displacement. The dominant elec- 
tron trajectory, shown as a dashed white line, crosses the 
lower edge of the Brillouin zone, but it doesn't reach its 
upper edge. Since Vc(k) = at the edges of the Brillouin 
zone, the electron is displaced most during the time when 
kx > 0. Consequently, the final electron displacement is 
positive, and the transferred charge is negative. 

Note that the temporal evolution of the current in 
Fig. [8] contains frequency components higher than those 
of the driving field. The same effect in a different param- 



Figure 8: The effect of dynamical Bloch oscillations in a sim- 
ulation with the following parameters: £0 = 10 GVm~^, 
V5CEP = 0, r = 2.3 fs, K = 0.1 fs-\ 7 = 0.3 fs^^ (a) 
Time evolution of the conduction-band electron population 
nc{kx,t) for a laser field with an amplitude sufficient for ac- 
celerating excited electrons to the edges of the Brillouin zone. 
The white dashed line represent the solution of Eq. ([6| with 
k(0) = 0. (b) The electric field of the laser pulse £x{t), the 
induced current density jx and the transferred charge Q{t). 



eter regime was reported to contribute to the generation 
of nonperturbative high-order harmonics in solids fl3|. 

Figs. [H] and [5] demonstrate that DBOs have a big im- 
pact on the CEP dependence of the transferred charge 
Q(</'CEp)- In Fig. [HI we investigate it more systemati- 
cally by plotting the normalized transferred charge 



Q(</5CEP,fo) 



Q('<gCEP,^o) 
(?^c}max(^o) 



(11) 



for different laser intensities. The £q dependent normal- 
ization factor 



max 

t 



/g2 "c(k, fo, t, ipc-E-p = 0) (fk 



(12) 

which is plotted in Fig. [9][c), denotes the maximum of 
the time dependent average population in the conduction 
band for (^cep = 0- 

One immediately observes that the extrema of 
QifCEp) shift as the laser intensity increases. By in- 
specting these simulations, we identified DBOs for peak 
laser fields above 7 GVm~^. Below this limit, (3(v?cep) 
is approximately proportional to cos((pcEp), and only the 
amplitude of Q{^pc'ep) increases with the peak laser in- 
tensity. For intensities high enough to induce DBOs, 
the extrema of Q(<y5cEp) change their positions. For 
Eq = \b GVm^^, maxima and minima exchange their 
places. 

To relate these observations to the semiclassical anal- 
ysis, we focus on the cosine pulse and plot the electron 
displacement at fmax = 15 fs as a function of the peak 
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electric field in Fig.|9l[a). We see that this analysis quali- 
tatively explains our numerical results: Q(v?cep = Ojfo) 
and s{to — 0,£o) have their extrema at approximately 
the same values of the peak electric field Eq. 




-0.2 0.2 0.6 1 -3 -2-10 1 2 3 0.1 0.2 0.3 0.4 

Semiclassical displacement CEP (radians) <np>n,a^ x 10^ 

(arb. units) 



of charge carriers in a wide gap bulk dielectric exposed to 
intense few-cycle laser pulses. Our most important find- 
ing is that whenever a laser field drives electrons close 
to or beyond the edges of the Brillouin zone, this has a 
significant impact on CEP sensitive measurements like 
those reported in ■ This effect survives even if we ac- 
count for dephasing and assume electron-phonon scatter- 
ing rates as high as 7 ~ 10^^ s~^, which is close to scatter- 
ing rates reported for Si02 [lOiHII- Thus, the detection 
of the total transferred charge can be used to measure 
signatures of dynamical Bloch oscillations in bulk solids. 
Furthermore, we show that this effect can be qualita- 
tively explained by evaluating the semiclassical electron 
displacement for dominant electron trajectories. 



Figure 9: (a) The semiclassical electron displacement (|10p 
for (fcEP = 0, to = and fmax = 15 fs. The parameters 
of the laser pulse are r = 2.3 fs and ujq = 2.51 fs~^. (b) 
The normalized transferred charge Q, defined by Eq. (|lip . 
as a function of of the CEP and the amphtude of the laser 
pulse. For dephasing and relaxation, we used k — 0.1 fs~^ 
and 7 = 0.3 fs~^. (c) The maximal average population in 
the conduction band (nc)max(fo), defined by Eq. (fT2)) . for 
V?CEP = 0. 



IV. SUMMARY 

Using an admittedly simple model, we have investi- 
gated the multiphoton injection and laser-driven motion 
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